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We apply the Callan-Harvey anomaly inflow mechanism to the study of QCD (chromoelec- 
tric) flux tubes, quark (pair)-creation and chiral magnetic effect, using new variables from the 
Cho-Faddeev-Niemi decomposition of the gauge potential. A phenomenological description 
of chromoelectric flux tubes is obtained by studying a gauged Nambu-Jona-Lasinio effective 
Lagrangian, derived from the original QCD Lagrangian. At the quantum level, quark con- 
densates in the QCD vacuum may form a vortex-like structure in a chromoelectric flux tube. 
Quark zero modes trapped in the vortex are chiral and lead to a two-dimensional gauge 
anomaly. To cancel it an effective Chern-Simons coupling is needed and hence a topological 
charge density term naturally appears. 
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I. INTRODUCTION 

Flux tubes in quantum chromodynamics (QCD) play important roles in many interesting 
places like confinement, quark pair creation, hadron structure and phase transitions. In the dual- 



superconductor description of QCD vacuum [l|, the chromoelectric field lines between color sources, 
like a quark and antiquark pair, are squeezed into a narrow flux tube along the line connecting the 
pair. The resulting potential is linear and consequently leads to color confinement. On the other 
hand, these flux tubes can provide chromoelectric field strong enough for quark pair creation, via 
the usual Schwinger mechanism (note that this is a non-perturbative particle creation mecha- 
nism) [s]. QCD fiux tubes are also well-studied in the formation of quark gluon plasma (QGP) in 
heavy-ion collisions. 

Another intensively studied object, is the topological charge and its distribution in the QCD 
vacuum. Different topological defects in the QCD vacuum have been explored by large-A'^c meth- 
ods, holographic QCD and lattice simulations. Among different topological defects vortices are 
particularly interesting, especially when we consider its possible connection the fiux tubes and how 
it interacts with the quarks. For example, the localization of quark zero modes on the topological 
defects leads to an interesting phenomenon — anomaly infiow jj] which could shed some light on 
the study of topological charge distribution in the QCD vacuum. It has been shown J, 5] that 
chiral fermion zero modes are localized on an axion string or a domain-wall embedded in a higher 
dimensional spacetime. When coupled to some external gauge potential, gauge anomaly appears 
on the string or domain-wall and it is cancelled by the gauge variation of an effective action, a 
Chern-Simons type coupling living in the bulk of the higher dimensional spacetime. 

The purpose of this paper is to provide a new description for the QCD fiux tubes and to study 
how topological charge emerges via the anomaly infiow mechanism. Our formulation is based on 
using a new set of variables for the non-Abelian gauge theory (Cho-Faddeev-Niemi) P, Q], i-e. the 
usual gauge potential A^^ splits into two new gauge potentials and B^^, which are convenient and 
efficient in studying both confinement and chiral symmetry breaking. As it has been shown in Ref. 
a, Isj, the gauge potential alone is responsible for the Wilson loop and the Polykov loop at the 
operator level, while the potential may be integrated out in the presence of dynamical quarks 
to produce a nonlocal four-fermion Nambu-Jona-Lasinio (NJL) interaction, which can be used 



to study chiral symmetry breaking as an effective Lagrangian method This decomposition 

and the resultant physics are more complicated in our formulation, since it describes the gauge 
configuration of a fiux tube with strong chromoelectric fields (plus chromomagnetic fields). In fact 
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if the flux-tube field lines are strong enough the chiral symmetry can be restored {l^ . Therefore the 
quark condensates inside the flux tube vanishes, i.e. {qq) = 0, while outside the flux tube the QCD 
vacuum has {qq) 7^ 0. Furthermore if this configuration carries a nonzero topological charge, it is 
possible for the quark condensates to form a vortex with a nonvanishing winding number, or even 
vortex lattice. We then can show that the chiral anomaly emerges in the bulk as a Chern-Simons 
effective action term. When chromomagnetic fields are included, this formulation can also be used 
to describe the chiral magnetic effect in QGP [15]. 

This paper is arranged as follows: Sec. II is to introduce the Cho-Faddeev-Niemi decomposition 
of non-Abelian gauge theories; Sec. Ill has a brief review of the Callan-Harvey anomaly inflow 
mechanism; We then study chromoelectric flux tubes in details in Sec IV; As an application quark 
creation due to chiral magnetic effect are considered in Sec. V; Discussions and conclusions are 
given in Sec. VI. 



II. CHO-FADDEEV-NIEMI DECOMPOSITION 



In this section we briefly review the reformulation of a non-Abelian gauge theory in terms of a 



set of new variables, or Cho-Faddeev-Niemi decomposition 6l4lCll|. We consider the SU(2) ca,se for 
simplicity and the generalization to SU(3) or SU(N) can be found in [3, S]- As shown in jol this 
decomposition is quite convenient in studying the crossover of the confinement and chiral symmetry 
creaking and also the topological feature of the gauge configurations. The decomposition reads 

mm 



(1) 



where the new variables and i?^ are defined as 



\ = {Af, ■ n)n + ig ^ [n, d^n] 
Bf, = zc/"^[V^n,n] 



(2) 



with V^n = — ig\A^^ n] and a unit vector n in the color space, satisfying 



n • n = 1 



(3) 



One can define a field strength 



(4) 
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in terms of A^^ and n 

F^,u = G^i^n, (5) 
Gf,u = d^{A,y ■ n) - d,y{Af, ■ n) + ig'^n ■ [d^n, d„n] 

The original field strength can be decomposed as well 

-^/tti' ~ -^jlu ~^ Ffiu (6) 

where fJj^ and are defined as 

-^i = (Ff^u ■ n)n, 

F^i = V^i?, - (7) 

with = d^ — ig[A^, ]. To see the Abelian character of this decomposition, we follow Ref. [3] 
and introduce a basis (n^,n^,n^) in the SU{2) color space. The color vector n can be identified 
with n'^. Then the gauge potential A^ can be rewritten as 

A, = n- + iC^ + Hf,)n' (8) 

where fi™ represents a classical QCD vacuum 

= -H^n'^ (9) 

since the corresponding field strength vanishes. The chromomagnetic potential is defined by 

H^^ = -^e'''''n%n^ (10) 

and the chromoelectric potential = A^^ ■ n^. The field strength F^^^, becomes 

Ff,, = {C^, + H^,)n^ (11) 

where the chromoelectric field strength C^u and the chromomagnetic field strength Hf^u are given 
by 

C^. = d^C, - d,C^, (12) 

H^, = d^H'i - d,Hl + ■ [9^, a,]n2 (13) 

respectively, where the last term vanishes except for the region where the topological defect is 
located, as it can be seen in the next section when we consider a stringy configuration. The action 
of pure Yang-Mills theory then becomes 



Sym = J d^'xl-^F^ 



4 fj-i^i 
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where B^i, = —ig[B^,B^] and we have used the reduction condition 



(14) 



V^B^' = 



(15) 



This condition also makes the dynamical degrees of freedom of the new variables match with those 
of original gauge potential. It is useful to rewrite the action in the form of jjo] 



Sym = J - -G^^ - - -B^'Q^^B'^ 



with the operator Qnu defined as 



(16) 



(17) 



The quark part of the QCD action in terms of the new variables will be given and studied in Sec. 
IV. 



III. ANOMALY INFLOW MECHANISM 



The localization of chiral fermion zero modes on vortices or domain-walls brings an interesting 



topic — anomaly inflow mechanism \^\ into the study of topological charge distribution in the QCD 
vacuum. It has been shown J] that chiral fermion zero modes are localized on an axion string or a 
domain-wall embedded in a higher dimensional spacetime. When coupled to some external gauge 
potential, gauge anomaly appears on the string or domain-wall and it is cancelled by the gauge 
variation of an effective action, a Chern-Simons type coupling living in the higher dimensional 
spacetime. In general, chiral anomaly in (2n -|- 2)-dimensions specifies the coupling between the 
gauge fields and the vortex. This coupling leads to the net flux of gauge charge into the vortex, 
cancelling the gauge anomaly due to the chiral zero-mode fermions in 2n-dimensions. This is how 
the chiral anomaly in (2n -|- 2)-dimensions is related to the gauge anomaly in 2n-dimensions. For 



geometric point of view, they are related by the Stora-Zumino descent equations 

Let us take Callan and Harvey's axion-string model as an example — an axion string, 
described by a complex scalar field $ = f{p)e^^, is embedded in a four-dimensional bulk space. 
Fermions coupled to the string have chiral zero modes localized on it. If the fermions are also 
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coupled to some gauge potential, there is a two-dimensional gauge anomaly due to the fermion 
zero modes 



1 



D''Jk = —e'^diA., i,j,A: = 0,l 



27r 



a 



which can only be cancelled by a Chern-Simons coupling similar to ()20p [4| 

1 



2 , (fxa^eKf" 



where is the phase angle of the axion field The effective coupling 



(18) 



(19) 



(20) 



contains the Chern-Simons current 



(21) 



which comes from the Chern-Simons form JCcs = Tr(A A F — ^A A A A A). Integrated by parts it 
gives a term resembling the QCD 0-term 



(22) 



with the help of the identity 



dJCr. = Ti F AF. 



(23) 



Note that the topological charge density (\22h is gauge invariant, while the current-current coupling 
(j20p is not. This is because the Chern-Simons form transforms under the gauge transformation 



A^A' = g ^Ag + g ^dg, 
JCcs Ks = ^cs - dTx{dgg-~^A) - h:T[{g-^ dgf]. 



(24) 



The gauge variation of (120 p is crucial in the anomaly-inflow mechanism, as it can be seen from the 



descent equation 



13|, 



5 [ d'dAK,l=! d^AdlCl 
Jm Jm 



/Co 



(25) 
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which relates the four dimensional axial anomaly dKS and the lower-dimensional gauge anomaly 



IC2 on the topo. 
cPi!} is singular 



ogical defect S. Note that the phase angle 'd is ambiguous at the origin and hence 



(fd = 27r 5 (x) 5 {y)dx A dy, (26) 
hence we cannot apply the identity = in (j25p . There is a subtlety in using covariant anomaly 



or consistent anomaly and it has been clarified in [5!]. We will show how this anomaly inflow 
mechanism on the string-like topological defects can be applied to the QCD flux tubes in the 
next section. Similar applications for the membrane-like topological defects can be found in Refs. 
Q) where we have studied the QCD membranes. 



IV. CHROMOELECTRIC FLUX TUBE AND ANOMALY INFLOW 

Chromoelectric flux tubes can be generated by quark pair creation, as described in Ref. [s^, in 
the process of hadron production in e~^e~ annihilation. Given enough energy, a quark pair qq is 
created at some point and then move away in opposite direction at a speed v close to the speed 
of light {v ~ c). The chromoelectric flux tube is then formed when the distance between the 
quark pair becomes 1 GeV~^ or more. New quark pairs can be produced in a cascade way. The 
confining constant chromoelectric field inside the flux tube pulls the new quark pairs apart until 
shorter flux tubes with lower energy form. 

We first consider a simplified case: There already exists a static pair of quark and antiquark 
qq and a chromoelectric flux tube connects the qq pair. We also assume that the flux tube is 
stable again the quantum fluctuations. These are possible scenarios: Firstly the chromoelectric 

flux tube is treated just like an ordinary electric flux tube, although the color degree of freedom 

n 

has to be included. This has been studied in details in |3] and other references, therefore we focus 
on another scenario, where the quark condensates are taken into account, and their interactions 
with quarks are studied as in the usual Nambu-Jona-Lasinio (NJL) effective methods in chiral 
symmetry breaking and other issues. Just like a magnetic monopole-antimonopole pair creates a 
string-like region of normal conduction in a superconductor, the quark-antiquark pair creates a 
region where the chiral quark condensate vanishes, i.e. (qq) = 0, while outside the chromoelectric 
flux tube the QCD vacuum has (qq) ^ 0. This is easy to understand since in sufficiently strong 
chromoelectric field, the qq pairs can be pulled apart and hence the quark condensates disappear. 
As it will be shown later in this section, it is possible for quark condensates to form a vortex in 
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the chromoelectric flux tubes. This is similar to the phenomenon in condensed matter physics that 
rotating Bose-Einstein condensates can form vortices when the rotating speed goes beyond some 
critical value. 

Using the new variables introduced in Sec. II, the quark part of the QCD action becomes 

5c,ua.k = J d^xli^ii^^^V^ - Mq)iP + 5^7^ia^S^] (27) 
Integrating out yields the gauged NJL effective action 

[ij{x + y)TA^{x - y) ^{x - y)rAHx + y)]) (28) 

where Ta are matrices in Dirac, color and flavor spaces. The non-local structure of the NJL 
coupling is not needed in this paper and we only consider the standard (gauged) NJL limiting case 
where G{x) = G6{x) and G is the usual NJL four-fermion coupling constant. The effective NJL 
action, which must keep the original symmetry of QCD action, can be parametrized as 

-Ceff = ^ [^7^(5^ - igAf,) - G{a + ifn) + • • • ] V (29) 

where the ellipsis includes contributions from the vector and pseudovector channels. So far this 
is quite similar to the usual NJL formulation except the "gauged" piece, i.e. the terms with the 
gauge potential Afj_. The usual NJL Lagrangian (perhaps even with some extensions) does not 
lead to a dual superconductor picture so we have to assume that the A^-related abelian dominance 
lead to quark confinement. As mentioned before the Wilson loop operator and the Polyakov loop 
operator only depend on but we will not consider the confinement problem in the present paper. 
Here we only need assume that the chromoelectric field lines are squeezed into tubes and hence 
the force between a quark-antiquark pair is linear. For the quark condensates to form a vortex, 
$ = cj + zvr = /(p)e*^ should satisfy a boundary condition 

/(O) = 0, /(oo) = constant (30) 

and this suggests that we choose 



|$|(0) = vV+^ (0) =0, 1^1 (oo) = constant (31) 
Can these conditions be satisfied? Follow the standard process of computing the effective potential 

n 

via Schwinger's proper time approach [2], one obtains 

Vefr(cT,^) = -^(^' + - ^Trln ({8^ - igA^f - ^a^,F^, + + vr^ - ie] (32) 
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For constant (chromo) electric and/or (chromo)magnetic background, this effective potential has 



been calculated before (see [12] for example) and one of the conclusions is that when the chro- 
moelectric field exceeds some critical strength, say Ecrit ~ 4GeV/fm [12'], the chiral symmetry is 
restored inside a chromoelectric flux tube, and we then have a normal phase with chiral symme- 
try, similar to the normal phase inside the magnetic flux tube of a superconductor. Therefore for 
strong chromoelectric flux tubes, the boundary condition (]3ip can be satisfied and it is possible 
that the quark condensates form a non-trivial topological defect — a vortex at the same region 
where the chromoelectric flux tube is located. Moreover, vortex lattice can be formed in the QCD 
vacuum just like those Abrikosov lattices in type-II superconductors. A couple of remarks are in 
order: First, the boundary condition cannot guarantee the existence of vortex solution, as a trivial 
configuration with a vanishing winding number is still a possible solution; Second, these vortices 
emerges at quantum level, and they are not solutions of the QCD vacuum at the classical level, 
which may not exist. 

Now let us study the effective Lagrangian (|29p in the background of ()3ip and a longitudinal 
gauge potential 

4 = (io,^,0,0), (33) 

for example, = (0, —Et, 0, 0) corresponds to an electrical field E background along the xi-axis 
(or the z-axis) and we will consider more general configurations in the next section. Writing ip in 
terms of ipR^L = 1/2 (1 ± r^) ip, the Dirac equation splits into two equations, similar to J] 

h\di - igAi)iPL + cos 9 + -/^ sme)dpipL 

+ fip)e-''^R = 0, 
if{di - igAi)^}^ + z(7^ cos + 7^ sin 0)dp'4}R 

+ f{p)e+''^L = 0, (34) 

which can be converted to equations of motion for a two-dimensional spinor Xl{xq,xi) 

i-f\di - igAi)xL = (35) 
and XL satisfies {—r^r^)xL = —XL- Then solution to ([31]) obtains the same exponential profile as 



in Ref. 



^L = XL exp [ - / f{p')dp'] 
Jo 



and ipR = i^'^ipL- Note that the difference between equations (p4]l . (|35]) and those in Ref. 



(36) 



i is 



that we have included the gauge potential A^ explicitly. The exponential profile in ([36|) shows 
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that chiral zero modes are localized on the flux tube. However, let us keep in mind that this is 
based on the assumption (I33p . Other gauge configurations could modify the exponential profile or 
even remove it, which indicates the transverse motion of quark zero modes. Here there is also a 
two-dimensional chirality is defined by the — 7*^7^ which is the two-dimensional analogue of 7^ in 
four dimensions. From the eqt. ()35p we see that chiral zero fermions are coupled to a gauge field 
and it is well known that a gauge anomaly appears in the flux tube 
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m 



V''Jk = ^e'^d,Aj, i,j,k = 0,1. (37) 

This seems to be a problem since the original theory QCD should not have gauge anomaly in the 
four dimensional spacetime. The gauge anomaly on the vortex configuration is actually cancelled 
by a topological term. 

The gauge potential Jl^ can be rewritten as [3] 

A^ = % + Z^n (38) 

where Z^j, = ■ n + is the U{1) potential introduced in ([8]) [ifl] and the corresponding field 
strength 

F^u = Z^^n (39) 

Z^, = d^Z, - d,Z^ (40) 



which shows the abelian character of Z^. As it is calculated in jj], the induced current in the QCD 
vacuum is 

1 



and it leads to an effective action 



167r2 

where the Chern-Simons current is 



■^m" = ^'>^^prZ'^9'^(^ (41) 



>Cc.e = -T-^ / d^xd^OK^' (42) 



= e^.prZ'ZP- (43) 



It is easy to see that the gauge variation of £(.-s cancels precisely the gauge anomaly on the vortex 
configuration. 
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Therefore we have shown that a strong chromoelectric flux tube may correspond to a vortex 
configuration of quark condensates. Quark zero modes located at the vortex becomes chiral which 
leads to gauge anomaly in the two-dimensional defect. This can only be cancelled by a topological 
term (Chern-Simons) in the bulk of the four-dimensional spacetime and hence relates the gauge 
anomaly in the vortex to the chiral anomaly in the bulk. From the current inflow point of view 
one may consider the flux tube as an "anomaly battery". It is charged by the current (]4ip whose 
only non-vanishing component is Jp = gE/{4:'iTp) (see Fig 1). It is also interesting that this current 
flows in the direction perpendicular to the chromoelectric field E like that of the Hall effect, as 



first noticed by Callan and Harvey 




FIG. 1. The anomaly- inflow picture of a chromoelectric flux tube in the QCD vacuum 



V. QUARK CREATION AND CHIRAL MAGNETIC EFFECT 

In this section we use our formulation to describe an interesting phenomenon in quark creation 
— the chiral magnet effect [l^. We will skip the chirality-conserving process since this is just the 
usual pair creation in strong (chromoelectric) fields via the Schwinger pair-creation mechanism 
Q], where the anomaly does not seem to play a role. The chiral magnet effect is different. The 
formation of QGP also involves color flux tubes. During the ultra-relativistic collision of heavy 
nuclei, net color charges develop and consequently color flux tubes connect these color charges. 
Quark-antiquark and gluon pairs are produced from the color fields via the Schwinger mechanism 
as mentioned before. Pair creation can happen in chromoelectric flux tubes and be enhanced when 
(chromo) magnetic field has a non-vanishing component in the same direction of the chromoelectric 
field, more precisely, E ■ B ^ 0. This is called chiral magnetic effect [3]. There are different 
approaches in understanding chiral anomalies, e.g. perturbation theory, path integral and index 
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theorem. The pair creation explanation in [3] is appeahng and particularly suitable in our case. 
From eqt. ([35]) and = {0, -E{t), 0,0) 

iY{d^-igZi)xL = (44) 

which can be readily converted to 



For a uniform chromoelectric field E(t) = Et, this reduces to the case studied in 



18( 1 one can 



apply a very simple argument which leads to the axial anomaly: the creation/annihilation rate 
of right/left-handed particles per unit time and unit length depends on the change of the fermi 
surface. According to {isl the annihilation rate of the left-handed particles is 

based on the dispersion relation from (I44p u}{p) = —p and the classical equation of a charged 
particle in a uniform electric field dp/dt = gE. Therefore eqt. (j46p gives the axial anomaly in 
two-dimensional spacetime 

nL-nR = £ f d^xeijZ'^ (47) 



and the axial anomaly is 



d'^jl = £e.,Z^^ (48) 



Notice that in the uniform (chromo)electric field case, 

^e,,Z'^ = E (49) 

This can be generalized to the four-dimensional spacetime. Assuming that at t = —oo the 
number of right-handed quarks is equal to the number of left-handed quarks. If the winding 
number of the vortex configuration is n, we then should have n chiral quarks created inside the 
flux tube at i = oo . Again it will lead to the gauge anomaly and need the anomaly inflow. This 
is more like a "recharging" process for "chiral batteries" , which was first discussed in [l^ . 

To study the chiral magnet effect we add a constant color magnetic field in the z-direction. 
Take the symmetric gauge = (0, —Et, —^By, ^Bx,0) and solve the Dirac equation in the new 
background with both constant color electric and magnetic fields 

il\di - igAi)iPL + i(7^ cos 9 + -/^ sin 9)dpipL 
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+ -gBp{j^ cos e - 7' sin 0) + fip)e-'^i^R = 0, 
iY{di - igAi)tl;R + i{-f^ cos O + j-^ sin 9)dpTpR 

+ ^5^P(7' cos 9 - 7' sin e)V^fi + f{p)e+''iljL = 0, (50) 

The profile of tlie quark zero modes now have an extra -B-dependent factor 

j' f{p')dp' - \gBp^] (51) 

( More general cases like time-dependent E(t),B{t) fields can be studied in a similar way). Now 
the direction of color magnetic field in the z-axis matters and it competes with the suppression 
factor due to the f{p) factor. Therefore we cannot treat the configuration as a string and must 
solve the equation in the four-dimensional spacetime. However, it is easy to see that if color electric 
and magnetic fields are (anti)parallel, the topological charge density 

^Zp.Z^^-' = B-E^O (52) 

The winding number of a gauge configuration 

Q = ^l d^xZ.^Z^^^ + (53) 

and hence the chirality is not conserved 

Nr-Nl = -2NjQ / (54) 

In order to have a complete analysis we need consider strong or moderate chromomagnetic fields 
and the initial states could be the lowest Landau level but this is beyond the scope of the present 
paper. 



VI. DISCUSSIONS AND CONCLUSIONS 

We close by pointing out that there are other formulation, e.g. the dual formulation in 

m 

= -C'/^n - g'^n x d^n + 0i d^n - (/<2 n x d^n (55) 

where (j)i^2 are two real scalar fields. The degree of freedom of the gauge potential can be moved 
around among the fields n and B^ and {C^,(j)i,(f)2) and n are dual to each other. With proper 
vevs for the above quantities, a (dual) abelian-Higgs model can be constructed [3]. It would be 
interesting to see how this is connected to our formulation. Note that these variables appear at 
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the classical level while some dynamical degrees of freedom in our model emerges at the quantum 
level. 

To summarize, we applied the anomaly-inflow mechanism to the studies of chromoelectric flux 
tubes, quark (pair)-creation and chiral magnetic effect, using new QCD variables from the Cho- 
Faddeev-Niemi decomposition. At quantum level, quark condensates in the QCD vacuum may 
form vortex-like structure at the region where a chromoelectric flux tube is located. Quark zero 
modes trapped in the vortex are chiral and generate two-dimensional gauge anomaly. To cancel 
this anomaly an effective Chern-Simons type action and hence a topological charge density term 
naturally appears. This scenario can also be used to study anomaly enhanced quark-pair creation, 
chiral magnetic effect and etc. In this paper we only considered straight string-like color flux 
tubes and their descendants. It is amusing to study more complicated cases, like the A-shape and 
Y-shape flux-tube configurations for baryons. 
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